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Abstract. Based on the CR formalism of algebraically special spacetimes by Hill, 
Lewandowski and Nurowski, we derive a nonlinear system of two real ODEs, of which 
the general solution determines a twisting type II (or more special) vacuum spacetime 
with two Killing vectors (commuting or not) and at most seven real parameters in 
addition to the cosmological constant A. To demonstrate a broad range of interesting 
spacetimes that these ODEs can capture, special solutions of various Petrov types are 
presented and described as they appear in this approach. They include Kerr-NUT, Kerr 
and Debney/Demiahski's type II, Lun's type II and III (subclasses of Held- Robinson) , 
MacCallum and Siklos' type HI (A < 0), and the type N solutions (A 7^ 0) we found 
in an earlier paper, along with a new class of type II solutions as a nontrivial limit 
of Kerr and Debney's type II solutions. Also, we discuss a situation in which the two 
ODEs can be reduced to one. However, constructing the general solution still remains 
an open problem. 

PACS numbers: 04.20.Jb, 02.40.Tt, 02.30.Hq 



Submitted to: Class. Quantum Grav. 



CR Structures and Twisting Vacuum Spacetimes: Type II and More Special 2 
1. Introduction 

All algebraically special Einstein spaces — vacuum but possibly with a nonzero 
cosmological constant A — possess a repeated principal null direction, which generates a 
foliation of the spacetime by a 3-parameter congruence of shearfree and null geodesies 
[H [2]. This 3- dimensional parameter space can be identified with a 3-dimensional 
real manifold described by the theory of CR structures with one complex and one 
real coordinate. CR structures were first introduced by Poincare and extensively 
studied by E. Cartan [3l H] . Good sources of background on the relationship between 
spacetimes and CR structures can be found, for instance, in the thesis of Nurowski 
[5] • Recently, Hill, Lewandowski and Nurowski [6] generalized earlier work of [TJ Ej to 
provide a new formulation of twisting algebraically special spacetimes with cosmological 
constant. It allows a classification of algebraically special spacetimes according to 
Cartan's classification of 3-dimensional CR structures. 

There is indeed a method to determine the equivalence of two solutions of Einstein's 
equations without having to construct explicit coordinate transformations that map one 
into the other, the idea of which was originated in work of Cartan and pushed forward by 
Brans p], Karlhede [9], and Skea [inilll]. Cartan also created a method for determining 
equivalence of two CR manifolds, which is much simpler than the method mentioned 
above for 4-dimensional manifolds. Because of the correspondence between the two, 
it allows a considerably simpler approach to determine the equivalence of two twisting 
algebraically special spacetimes, or, more usefully in this paper, the lack of such an 
equivalence, thereby guaranteeing that two solutions are distinct. In addition to this 
important reason for using CR structures, it also provides a different formulation of 
Einstein's equations, which exhibits certain invariant features that are desirable for 
calculations, as compared to other formulations. It is favorable to have this invariant 
approach to study the decomposition of Einstein's equations into some manageable 
form, and to have one that prefers non-zero values of the twist, since only a very limited 
number of such spacetimes with non-zero twists, in Petrov types II, III, and N, are 
actually available for study. 

In the theory of exact solutions, Einstein's equations are usually solved under 
assumptions of the existence of some symmetry group [12], i.e.. Killing vectors. 
For instance, Kerr and Debney [13] have determined all diverging (twisting or not) 
algebraically special vacuums (A = 0) with three or more Killing vectors. The case with 
two Killing vectors, however, is still not solved completely. We intend to address this 
problem in this paper with the extension to include a nonzero cosmological constant. 

Building on the work of Hill, Lewandowski, and Nurowski [6], we first present, in 
Section 2, the twisting type II vacuum metric formulated according to CR geometry, 
together with our calculated Weyl scalars. Then in Section 3, we establish the 
transformation from the CR formalism to the canonical frame that is widely used, 
e.g., in [12]. In Section 4, we generalize the ansatz that was found for twisting type N 
solutions in [13], thereby reducing the field equations to only two coupled real ODEs for 



CR Structures and Twisting Vacuum Spacetimes: Type II and More Special 



3 



two unknown functions of a single variable. From Section 5 to 8, we show, by using the 
transformation in Section 3, that a large variety of previously-known, twisting solutions 
of types II, III, as well as D ([I2], Chapters 29 and 38), with at least two Killing vectors, 
correspond to special solutions of these ODEs. Moreover, we study a special case when 
the number of ODEs can be reduced to one, which generalizes our previous results on 
type N. Though the general solution is yet to be found, we believe that this extension 
of ODEs is quite worthwhile and should provide a start for future steps forward in the 
study of twisting exact solutions. 

2. CR structures and the field equations 

A CR structur^ is a 3-dimensional real manifold M equipped with an equivalence class 
of pairs of 1-forms A (real) and /i (complex) satisfying 

A A /i A 7^ 0. 

Another pair (A',/i') is equivalent to (A,yu), iff there exist functions / 7^ (real) and 
h ^ 0, g (complex) on M such that 

A' = /A, fi' = hfi + gX, Ji =hji + g\. 

For our purpose, we further assume [6l [3] 

^ = dc, ji = dC, (1) 
dA = i/i A /i + (c/i + c/i) A A, (2) 

where C, and c are some complex- valued functions on M . Taking the closure of ([2]), we 
obtain a reality condition on the derivatives of c: 

dc = do, (3) 

The same function c also appears in the commutation relations of the dual basis of 
vector fields: 

{do,d,d) dual to (A,/i,/i), 

[d,d] = -ido, [do, d] = cdo, [do,d] = cdo. 

There is then the following theorem telling us how to construct an algebraically special 
spacetime on the basis of M. 

Theorem 1. |6l [15] The CR structure ([T]|l]) on M can be lifted to a spacetime 
Ai = M X M. equipped with the metric 

g = 2 {e^e'^ + e^e'^) , e^ = Vfi = P, 

e^ = VX, e^ = V (dr + W/i + W/i + ?/A) , 

where V ^ 0, H (real) and W (complex) are arbitrary functions on Ai. The spacetime 
([5]) admits a geodesic, shearfree and twisting null congruence along the vector field dr 
(r G M), of which the 3-parameter leaf spaces (r = const.) have the same CR structure 



^AU our considerations are local. 
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as M. It further satisfies the Einstein equation Ric(g) = Ag, iff the metric components 
can be written as 

P 



V 



>V = ia(e-"' + ll 



cos 



n 



H = -e2'^ + -e-"'^ + ge^^ + ge- 
a = c + 291ogp, 



n 



-2ir 



(6) 

(7) 
(8) 

3n + n 2 2 , '^dpdp-p [ddp + ddp) i - 

+ -Ap^ H ^ - - do iogp - do, (9) 



2p^ 



(10) 

- p^ 

where c, n (complex) and p (real), all functions on M (independent of r), satisfy the 
following set of equations: 



dc = do, 

- 1 3 
dd + dd + cd + cd + -cc + - (dc + dc) 

dn + Sen = 0, 
R33 = 0. 



n + n 2.0 
P = — 5- + ^Ap^ 



(11) 
(12) 

(13) 
(14) 



Here the Ricci tensor component R^^ as well as the Weyl scalars "^2, ^3 and \E'4 are 
given by 

R33= l^{d + 2c)[p'{dI-2A{2d log p + c)p') ] 



16 . 
+ — A 

P 



dd + dd + cd + cd + -cc + ^{dc + dc) ]p — _ -\p^ 



1 3 

2- i' 



po 



16i ^ / n 

\p^ 

n 



cos 



2p6 



e"- + l 



(15) 
(16) 



1^ [9J-2A(2aiogp + c)/] 



(17) 



+ 6i (29 log p + c) ^ (e^^^ - 1) - 4i9(^) (e''- + l)|e 



i^/^cos^f- 
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^4 = 2e-^"/2 cos^ (^0 I - (29 log p + c) [dl - 2A(2aiogp + c)p^] ^ ~ 

+ {3 + 2c) [dl - 2A(2aiogp + c)p2] + —dj 

2 - - - - n 

+-A \(d + c)(2d\ogp + c) + 2(2d\ogp + cf] - 3(2d\ogp + cf— e'^'' 

— — — ?2 — — / n 

(d — 29 log p) (29 log p + c) -TT + 3(29 logp + c)9( — r 

pb \pb 

n (169 logp + c)9n — 9^n 



+ 



3 (29^ logp — 16(9 logp)^ + 2c91ogp + c^) 



p6 p6 



+ 3 (29^ log^, - 8(9 log^,)^ + 8c91og^. - 9c - c^) ^ + n2mogp -c)9n - 29 n 

+2 (9^ logp - 2(9 log p)^ + 5c91ogp - 9c - 3c^) ^ + i^d log p - 5c)dn - d^n \ 
\ ' p° p'^ ) 

with the function I defined by 

/ = 9 (91ogp + c) + (91ogp + cf . 

Following the same procedure as [6], which uses Cartan's structure equations to 
calculate the curvature tensor, we present our calculated -R33, \E'3 and '^^ above with 
nonzero A and n, as a complement to [6]. Moreover, to facilitate future calculations, we 
have arranged the expression of -R33 so that its second square bracket can be immediately 
removed by the field equation (fT2!) . whereas the terms dl — 2A(291ogp + c)p^ are made 
prominent as they also appear in \&3 and \E'4. 

To solve the field equations (1111114^ in practice, one needs to introduce a real 
coordinate system (x, y, u) on M such that 

C = x + iy, = I - ' x ^ dw + LdC + i^dC , . 

9 = 9c-L9„, 9o = i(9L - 9L)9„, i{dL - dL) ' ^ ^ 

with L = L{(X,u) a complex-valued function [16] satisfying 

BL-dL^ 0, (19) 

which is needed for a nonzero twist (cf. (I27ll ). In addition, L relates to the function c 

by 

c = -9 ln(9L - dl) - 9„L, (20) 

as imposed by the commutation relations (j4]). Hence generally, the system ( fTTIfT4|) are 
in fact PDEs for the unknown functions L, n and p of the coordinate variables (C, C '^)- 
For other possible coordinate choices, the metric (fTl l20l) admits the following 
coordinate freedom ([IT], see Section 2.6): 

r' = r, C' = /(C), «' = F(C,C,«), d^F^Q, (21) 
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with f{() holomorphic and F a real- valued function, which generates the transformation 
laws 

/i' = /V, A' = /7'A, /' = d//dC, (22) 
5' = -^5, ^o = Jj;9o, [d',d']=-idi [d'„d'] = c'di (23) 

p-^p, n'^^^n, (24) 

L' = -jid^F ~ Ld^F) = -j^dF, B'L' -d'L' = ^{dL-dL). (25) 

Note that the function F does not appear above on the level of c, p and n, as well 
as d and 9o, the fact of which indicates an invariant feature of the CR formalism. As 
expected, the field equations for the new p', d and n' take on the same form of (11 111 141) 
with {do,d,d) simply replaced by {dQ,d',d'). These transformation properties will be 
used to simplify our metrics. 

3. Transformations to the canonical frame 

Given the algebraically special twisting metric form (fTl l20|) formulated according to CR 
geometry, it is important to know how it is related to other pre-existing formalisms that 
have been extensively studied in the past. Here we quote from [12] (p. 439-441) a most 
commonly used one by Kerr, Debney et al. [HI [191 [20] . Foi' simplicity, we only consider 
A = and follow closely the notation of [12] with sub- or superscript s added to avoid 
confusion. 

Theorem 2. A spacetime admits a geodesic, shearfree and twisting null congruence 
along the vector field dr^ and satisfies the Einstein equation Ric(g) = 0, iff the metric 
can be written as 

g = 2{ujW - a;V), uj^ = --^ = u\ 

PsPs (26) 

= du + LdC + LdC, oj^ = dr, + W^dC + W^dC + HsCO^, 
with metric components 

2iS 

= -(r, + iS,), -^ = BL-dLy^O, 
Ws = p-'duL + id^s. d = d^-Ldu, 
Hs = -Ks - r^du \ogPs ^2j^Y? ' 



= 2p2Re[a {d\ogPs-duL)] 
such that the functions m^, M^, (real) and L (complex), all only dependent on the 
coordinates (CjC"^); satisfy a system of PDEs: 

P:^Ms = Im ddddVs, Ps = duVs, (28) 
d{ms + iMs) = 3(m, + iM,)9„L, (29) 
du [P^'inis + iMs)] =Ps[d + 2(9 log - duL)]dIs, (30) 



CR Structures and Twisting Vacuum Spacetimes: Type II and More Special 7 



where the function Ig is defined by 

I, = 8 (aiogP, - d^L) + {BiogPs - d^lf = P;^duddVs. (3i) 

Additionally, the Weyl scalar \E'| is given by 
% = (m, + iMs)pl. 

In this metric form, the coordinates (CC"") cind the function L have been chosen 
identically with those introduced in (fT8l) : hence, each is not given a sub- or superscript 
s. Taking A = and by a tedious but straightforward calculation, one can show that 
the metrics ( 12611311) and (ITI l20|) are equivalent to each other by the transformation [17] 

^ \{dL - dl) ' ^^^^ 

= — 71^ — ^ tan f-^ , Irl < vr, (33) 
i{dL-dL) V2/ ' ' ' ' ^ ' 

16(n — n) 16(n + n) , 



with the inverse 



p='-{dL~dl)Ps, (35) 

r = 2arctan| — — ^^^1, (36) 

V i{dL - dL)P^ 7 ' ^ ' 

n=^{m, + iM,){dL-dLf. (37) 



In particular, the field equation ( 1281) can be transformed into ( U21) with A = 0, despite 
their drastically different appearances. Also, one gets = / when substituting fl32|) into 
the definition (13T|) . For more details about these transformations, one may see Section 
2.6 of [IT]. The relation ( 135|) and (1371) together with (120!) will be used later to translate 
known solutions of the canonical field equations (I28tl30p to solutions of ( ITTIfT^ . 



4. Reductions to ODEs 



Now we go back to the metric (fTl l20|) . Following the same idea as [T3| for solving the 
field equations, we assume that the unknowns p, c and n have no u-dependence, i.e., 
SqP = Oqc = Oqu = 0. This assumption avoids the involvement of the function L inside 
the operator d since now we have, e.g., dp = d(^p. Therefore the system (IlltiMp becomes 
effectively PDEs for the unknowns c, p and n instead of L, p and n. Once c = c(C, C) 
is solved, one may further determine a function L without w-dependence from (120|) (cf. 
(H9l) ). Altogether, this means that the resulting spacetime shall possess a Killing vector 
in the u-direction, which is, in fact, an assumption widely used in many research articles 
on algebraically special solutions (see, e.g., [12] Chapter 29). 

We apply the assumption and rewrite the system ( ITTIfMl) (likewise for and ^^4) 

as 

d^c = %, (38) 
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1 3 n + n 2 
2di;d^p + cdc^p + cd^p + -ccp + - {dc_c)p = — h -Ap^, (39) 

9^n + 3cn = 0, (40) 
((9^ + 2c) [/9^J - 2A(29f logp + c)p^] = 0, (41) 

where in the last equation {R^s = 0) we have used (1121) to simplify the expression of 
-R33, and the function / is given by 

I = {d^ \ogp + c) + (9^ logp + c)^ . 

This is the set of PDEs we aim to solve for the unknowns p(C, C); c(C, C) ^'^^ "^{Ci 0- 

Generalizing the ansatz [T^ [T7] we found from the classical symmetries [21] of the 
type N case of ( 138MT1) with = 0, we assume the following forms for the unknowns: 



p(C,C)^^, c,C,C)^M±M£>±^, ni,,0.?M±^ (42, 
with a new real variable 

idC - / =1^1 -AC- 



A ^ J A J A 

Here the function A = A{() is an arbitrary function of ( that is sufficiently smooth, 
and the constant Ci and the undetermined functions Fi_4{z) are all real- valued. The 
constraint equation (l38l) has been taken into account in the form of c(C,C) so that it is 
satisfied. 

Inserting the ansatz ( l42|) into ( !39ll4Tl) , we obtain a remarkable reduction to a system 
of four compatible real ODEs for Fi-4{z) only, with all other dependence on A, A 
factored out: 

= -F[' + F,F[ + ^AFf + i(3F^ - F| - C',)F^ + ^, (43) 

= F^- 3(^2^3 + C1F4), (44) 
= Fi + 3(CiF3-F2F4), (45) 
= {H' - 2F2H)' - 2F2{H' - 2F2H) + AClH, (46) 

where the function H{z) is defined by 

H = F{'Fi - {F[f - AFf - Fl^Fl 



The two inner equations are derived from the real and imaginary parts of ( l40l) . 
respectively. Note that they are linear in F3 and F4. Thus if F2 is given, one can 
solve them by 



F3 = exp j F2dz^ [Bi sin(3Ci^) - B2 cos(3Ci2)], 
F4 = exp j F2dz^ [B2 sin(3Ci^) + Bi cos(3Ciz)], 



(47) 
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where B12 are real constants. With F3 expressed in terms of F2, we are left with only 
two nonlinear equations for Fi and F2: 

= - Fi' + F2F[ + ^-XFl + i(3F^ - F| - Cl)F, 

+ M3C,z) - B2 cos(3C,z)], (48) 

0={H' - 2F2H)' - 2F2{H' - 2F2H) + AClH, 

which can be easily converted to a set of ODEs if one introduces, e.g., F2 = K'{z) to 
remove the integral in the first equation. The system psj) with f H7|) . or alternatively 
( 1431H6|) . constitutes the main result of this paper. 

Before moving on to solve ( 143H46|) in the next few sections, we should make a few 
general remarks concerning the metric (1T H2U]) equipped with the ansatz ( I42p . First, 
despite the appearance of a free function A{C,) 7^ in the ansatz, its different choices do 
not generate new metrics. In fact, using the coordinate change C = f '^q^^ permitted 
by (I2T]) and, accordingly, the transformation law with /' = 2/A{Q, one can always 
replace a function A by a constant A = 2. This is also consistent with the fact that the 
local CR structure determined by the function c in ( H2l) is independent of the choice 
of A{() [H]. Hence for simplicity, we can just set A{() = A{() = 2 without loss of 
generality, and hence z = lm( = y. With this choice of A, the function L and the 
1-form A can be determined from (!20l) and (ITSl) as 

L = -e-^''-" J ^^v(^J F^dz^ dz, (49) 

^ e^i^dw - 2 [/ exp (/ Fsd^) d^] dx ^^^^ 
exp (/ F2dz) 



Once again, though there exist other L's satisfying fl20|) . one can always use the 
remaining coordinate freedom u' = F{(,(,u) (cf. fl^T]) ) to convert them to the u- 
independent real expression (149|) [14j . Given such L and A, the class of metrics 
determined by (14211461) admits at least two Killing vectors 

Xi = du, X2 = d^- Ciu du, (51) 

with the commutation relation 

[Xi,X2] = —CiXi. 

These vectors, verifiable by direct calculation, are both inherited from the symmetries 
of the underlying CR structures [22|, [17] . 



5. Type D solutions: Kerr-NUT 

All type D vacuum solutions, twisting or not, are known [23l [2l], which include 
perhaps the most famous algebraically special solutions such as Kerr's rotating black- 
hole solution. Hence it is worthwhile to consider whether our equations (1431146 p can 
capture some of these physically important solutions. But first we should comment 
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that the CR formahsm (ITI I2D]) used here is constructed on just one single shearfree null 
congruence aligned with a multiple principal null direction , which is unique in type II, 
III and N spacetimes. However, a type D spacetime possesses two such congruences, 
each along one of the two doubly degenerate principal null directions, and consequently 
one cannot treat them at the same time in the CR formalism. This suggests that 
the CR formalism may not provide the most convenient approach for finding type D 
solutions, as compared to other approaches that are specially designed to make use of 
both congruences. 

A spacetime is of type D iff it satisfies the conditions 

3^3^4 - 2*^ = 0, *2 ^ 0. (52) 

This equality with the ansatz (H2l) applied gives rise to a number of lengthy ODEs (as 
one may sense by looking at the expression of \E'4) from the coefficients of various powers 
of e"' required to vanish. This fairly complicated situation (except when \l/3 = \l/4 = 0, 
see [Appendix A]) needs a specialized paper to elaborate; hence it is not further discussed 



here (also because there is no new type D solution to be found). Instead, by applying the 
results of Section [3l we will simply show that the Kerr-NUT solution can be retrieved 
as a special solution of (14311461) through the ansatz P2|) . 

In the canonical frame ( l26ll3T|) . the Kerr-NUT solution [12] (p. 453) is given by 

CC 

2 ' - (P^ 
m, = m, M, = M, A = 0, 



Ps = l + ^, L = -77^[2M + (M + a)CC], 



where m, M and a, each a real constant, are called the mass, the NUT parameter and 
the Kerr parameter, respectively. Then inserting them into fl5Sl) . fl2U|) and fl57|) . we 
obtain the following solution to (13 81141 P : 

4(M-a) + 2(M + a)CC 



P 



n 



(2 + CC)^ 
2C[2(M - 2a) + (M + a)CC] 
[2(M-a) + (M + a)CC](2 + CC)' 
16i[2(M - a) + (M + a)CC]^("^ + iM) 



(2 + CC)^ 

which, as expected, satisfies the condition (152|) for type D with nonzero \I'3 and 
Without a dependence on u, these expressions can be cast into the form of our ansatz 

m by 

A(c) = -iC, ^ = iog(CC), C'i = o, 

such that 

2[2(M-a) + (M + a)e^]e^/2 ^^{^m - a) + ^ae' - {M + a)e^^ 



(2 + e-)2 ' [2(M-a) + (M + a)e-](2 + e-)' 

16M[2(M-a) + (M + a)e^]V^ 16m[2(M - a) + (M + a)e^] ^ 

(2 + e-)9 ' ^~ (2 + e^)9 ' 
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One can verify that they are indeed a solution to (14311461) with A = 0. Note that the 
function A{() above will still be serving as a free function, as long as -Fi_4 are obtained. 
The Kerr-NUT solution is contained in the Demiahski solution as a special case (cf. 
(168|1 ). For other examples of type D solutions, see [Appendix B| and ( 169|) with B2 = 0. 



(54) 



6. Type N solutions with nonzero cosmological constant 

The type N solutions require 

^2 = ^3 = 0, ^4 ^ 0. 
These conditions lead to the following special case of fl43tl46p with F3 = F4 = 0: 

= -F^ + F^Fi + lAF^ + \{3F^ - F| - C',)Fr, 

^3(xH' - 2F2H - 2iCiH = 0, 
plus one inequality 

4A r 

^4 = - [2F^F[' + 6{F[f - 10(F2 + iC,)F,F[ 

- {F^ - 3F| - 6iCiF2 + 3C^)F^'j e"'"/' cos2(0 ^ 0' 

which imposes A ^ for type N. To better see that the system (ISlll is included in 
(14311461) . one can rewrite equation fH6|) as 

= {H' - 2F2H - 2iCiH)' - 2(^2 - iCi){H' - 2F2H - 2iCiH). 

In fact, the field equation R^-^ = can always be removed by \E'2 = = for general 
type N vacuums (cf. (I15m7p ). Notice that the second equation of (154|) is complex; we 
have two cases for solutions. 

Case 1: H = 0, A 0. This simpler case has been investigated in [13] (see [17] for 
more details). The equations for this case read 



= -F^ + F2F[ + lAF^ + \{3F^ - F| - Cl)Fi, 
= H = -F['Fi + {F[y + AF^ + F^F^. 
By introducing a real function J = J{z) and 



(55) 



Fi = ±y7, F2 = ^-AJ, J'>0, (56) 



2J' 

we can reduce the first equation of fl55|) to 

( T"V in 
J'" = ^±J- - 2AJJ" - —AiJ'Y - 2(AV2 + C^)J', (57) 
2(7 3 

while the second equation is automatically satisfied. Since this ODE has no explicit 
dependence on the variable z, we can immediately lower its order by the transformation 
J' = P{J)>0 such that 

A,o, (58) 
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which, in the case of Ci = 0, can be further reduced to an Abel ODE [23] of the first 
kind 

by J = exp(/ f{t)dt)/A, P{J) = texp(2 / /(t)dt)/A. Unfortunately, this Abel ODE has 
not been identified as a known solvable type. 

Various aspects of the equation fl58l) were examined in [T3|, including the weak 
Painleve property [26] and constructions of various special and series solutions. All 
degenerate solutions of type O (\l/4 = 0, conformally flat) were found. The only known 
type N solution with A 7^ in closed forms was first discovered by Leroy [27] and 
presented in the CR formalism by Nurowski [15]. It corresponds to 

P(j) = _iAj2_^>0, 
^ ^ 3 4A ' 

' 2ssin(|Ci(^ + Co))' ' tan(iCi(^ + Co))' 

with A = — < and Cq a real constant (removable by a translation z + Cq — )■ z). 
Nonetheless, the equation fl58|l does also admit type N solutions with A > and two 
additional parameters besides A and Ci. 

Case 2: if ^ 0, = 0, A ^ 0. The system ([54D becomes 

= -F'l + F^F[ + \KFl + 1(3F^ - 
= H' - 2F2H. 

Very little is known about the solutions of this system except for one of type O [17] 
given by 

^' = ^2s{flcoy ^' = ~VTc-o' ^ = As^{zlc,r ^^^^ 

with A = — < and Cq a real constant. Particularly, this solution has the 
hyperquadric CR structure (the most symmetrical one) [2B] . 

7. Type III solutions 

Similar to the case of type N, the equations for type III (\E'2 = 0, \E'3 7^ 0; ^3 = ^4 = 0) 
are given by 



= -F^ + F2F[ + lAF^ + \{3F^ - F| - C^)Fi, 
= {H' - 2F2H)' - 2F2{H' - 2F2H) + AClH, 



(63) 



which are subject to 

^3 oc if' - 2F2H - 2iCiH ^ 0. 

Using the first equation of ( l63i) . we can lower the order of the second ODE, such that 
the resulting set of equations contains derivatives up to the second-order in Fi and 
the third-order in F2. Hence the general solution carries another five real parameters in 
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addition to A and Ci. Considering that one of these parameters is simply the translation 
z ^ z + Cq (no explicit dependence on z in fl43ti46p ). and thus removable, one can see 
that the final type III metric determined by ( l63l) has at most six parameters including 
A and Ci (see the conclusions). 

Two classes of twisting type III vacuum solutions are known, respectively for A = 
and A < 0. The one with A = is due to Held [29] and Robinson [30], which generalizes 
the non-twisting Robinsion-Trautman type III vacuum solution and generally admits 
only one Killing vector du- The subclasses with two Killing vectors (commuting or not) 
were found by Lun [31]. It can be shown that Lun's case I type III metric corresponds 
to the following solution of ( 1631) : 



A = 0, Ci = 0, 



z-^ 

(5 - V^)Eiz2+v^/2 + (5 + v^)E2z2-v^/2 



Likewise, his case II type III metric puts forth a second solution of fl63l) : 
A = 0, Ci = i, 

_i 

-^1 = 77t{FiG + E2G ), 

4cos(f) 

VT3{E,G-E2G~') ,5^ (65) 
^ — ' tanf^ 



G{z) 



4cos(f)(EiG + E2G'-i) 4 

sm(|) + 1^^ ' 
cos(|) 



'.2n 



In both cases, £'1,2 are real constants. More details will be given in the next section as 
degenerate cases of the related type II solutions (1671) and ( ITOl) . 

The other known type III solution requires A < and is due to MacCallum and 
Siklos [32] (see also [12] p. 201). As a solution of ( j63l) . it is given by 



4:s{z + Co) 2{z + Co) 

with A = — s^, a real constant Gq and 

585 



*3 OC if ' - 2F2H 



256A(z + GoY 
8. Type II solutions 



Based on the structure of (14311461) . we can consider the type II solution (\E'2 7^ 0, 
3\E'2^4 — 2\E'| 7^ 0) according to three different cases. 
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Case 1: F3 = 0, F4 7^ ^ Ci = 0. The equations (1431146 p are reduced to 

= -F^ + F2FI + lAF^ + \{3F^ - F^)F,, 
= {H'- 2F2H)' - 2F2{H' - 2F2H), 

F4 = Biexpis F2dz 



14 



with Bi ^ a. real constant. Since the first two equations above are identical to (!63|l 
with Ci = (also cf. (155 p and (I6ip ). one can generate this kind of type II solutions 
directly from existing type N and III solutions, i.e., fl60|) with Ci = 0, f l6^ and f l66|) . 
which works as if one is adding a "mass source" to them. A similar idea can be found in 



12] p. 447. In addition, for type D solutions with F3 = and Ci = 0, see Appendix B 
Case 2: F3 7^ 0, Ci = 0. The associated equations are given by 







F3 



-F'{ + F2F[ + \KFl + i(3F^ - + -i 



= (i/' - 2F2H)' - 2F2{H' - 2F2H), 
F2 = 7^, F4 = Bi exp (3 / F2dz 



B1F3 



with Bi a real constant. Certainly one may use the third equation above to turn the 
first two into ODEs for Fi and -F3 only. 

Lun's case I solution [3T1 [33] with four parameters can be shown to belong to this 
case. It reads, in the canonical frame fl26]|3TD. 



' (C + C,f'\ 



T 3i _ 

L = X 

16 



( = x + iy, 

■'"^ [(3 + ^/U)Eix^'/^ + (3 - 
m,=m, Ms = M, A = 0, 
or, in the CR formalism as a solution of ( l38ll4T]) or ( H3M61 
A{C) = -2i, z = x, Ci = 0, 



13)E2X' 



13/2 



3iM 
32x^' 



P 



-,F.iz) - 
-lF2i^) 

{zG)~^ [(5 
1 



V3_ 

16^3 



G{z) 



13)^1^2+^/2 



27i 



n 



64 (^3W+ii^4(^))- 220,27/2 



- (5 + v^)E2^2-v^/2 + QMz^/^ 
G'\m + iM), 



(67) 



G{z) = 4 (Eiz2+VT3/2 + E2Z' 



13/2 



3Mz'^\ 



with m, M and Ei^2 real constants. When m = M = 0, the solution degenerates to the 
type III solution (16^ . 

The Kerr and Debney/Demiahski's four-parameter solution [12] (see p. 449) also 
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falls under this case. It is given by 

= 1 + y , L = - 



-iP 



A 



2M/C+(M + a)C + i6Clog 
0, 



nis = m, 
and corresponds to 

A(C) = -iC, ^ = iog(CC), C'i = o, 

[16(M - a) + 8(M + ay + 6Gi]e"/2 



4(2 + e^)2 

8[4(M - g) + 8ae" - (M + a)e^^] + 6^2 
8[2(Af - a) + (M + a)e^](2 + e^) + fed (2 + 
M[16(M - a) + 8(M + a)e^ + 

32(2 + e^)9 ' 
m[16(M - a) + 8(M + a)e-' + 



(6^ 



32(2 + e^)9 
Gi{z) = (2 + ln2 - z){2 - e") - 8, 
G2{z) = (3 + ln2-2)(4 + e2^)-8(ln2-2)e"-24, 

which is a solution of (1431146 p . Here m, M, a and b are four real parameters. Clearly, 
the Kerr-NUT solution ( j53l) is a special case with 6 = 0. 

Besides these known solutions, we have obtained an additional one (see the 
derivation in Case 3 below) which turns out to be a limiting case (Ci — > 0) of the 
solution ([72]): 

Ci = 0, A = 0, 

4^22 - C2 



C2Z + C3, F2 



(69) 



252^2 -C2^-C3' 

F3 = —B2F'l, F4 = BiF^. 

Here C2,3 and B12 are real constants. Its comparisons with Lun's and Demiahski's 
solutions will be discussed in Appendix C Particularly when B2 = 0, the solution 
becomes type D. 

Case 3: F3 7^ 0, Ci 7^ ^ F4 7^ 0. This corresponds to the most general case for 
solutions. As one may check, Lun's case II four-parameter solution [311 [33], which is 
given by 

Ps = \[liC + C)'/', C = x + iy, w = y/x, 

L = x''/^[l[E^iw + {W^ + l)V2)v^/2(y, _ v|3(^2 ^ 1)1/2) 

+E2{w + {w' + iy/'r^/\w + ^{w^ + 1)1/2)" 

+^[{m + iM)(l + iwf/^{2 - 3iw) + (m - iM)(l - iwf/^{2 + 3iw;)]|, 



— iMo 



(m 



+ iM)x3/2(l 



3/2 



A = 0, 
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can be converted to a solution of 



(70) 



^(C) = C, C = IC|e-^ Ci = l 

VSjE^G + E^G-^) 3V3 sin(f) (Msin(f) - mcos(|)) 
'~ 4cos(f) ^ 16cos5/2(|) 

_ 4cos3/2(|)[5sin(f )(EiG + EaG^^) + VT^iE^G - E^G-^)] - 3(MGi - niG^) 
^ ~ 16cos5/2(|)(EiG + E2G-i) +6sin(z) (Msin(f ) - mcos(f )) ' 

_ sV^Ff (Mcos(f ) +msin(f )) _ sV^Ff (-Msin(f ) + mcos(f )) 

" 28 C0S9/2(|) ' " 28 COs9/2(f) ' 

G{z) = [-i^Ey) ' G^i(^) = iMt) - |sin(f ) -5sin(|), 
G2(^) = |cos(f) + |cos(f)-5cos(f), 

with -Ei^2; M and m real constants. In the case of M = m = 0, the solution reduces to 
the type III solution ( l65|) . 

Besides Lun's example, we have also considered the special case of if = for the 
system (14311461) . which turns out to be fully soluble when A = 0. The derivation follows 
closely the type N Case 1 (cf. (l55l) ). and utilizes the same ansatz ( l56l) that makes the 
function H{z) vanish (hence, ( H6|) satisfied). More specifically, we have 

7" F' 
F2 = — - A J = ^ - AJ, J' > 0, 
2J' Fi 

F3 = F^exp (^-3A j Jdz^ [Bi sin(3Ciz) - ^2 cos(3Ci^)], (71) 

F4 = Fj^exp (^-3A j Jd2^ [B2 sin(3Ci2) + Bi cos(3Ci2)], 

with the last two equations derived from P7|) . Therefore when A vanishes, the equation 
( 143|) can be reduced to a linear ODE for Fi{z) alone: 

= -GlFi + 4[5i sin(3Ciz) - B2 cos(3Ci2)], A = 0, 

which has the general solution 

Fi = C2 cos(Ci(2 + Co)) - ^[fii sin(3Ci2) - B^ cos(3Ci^)], Ci ^ 0, (72) 



Fi 


= ±\ 




F3 


= F^ 


exp 


F4 


= F^ 


exp 



or, if Gi vanishes (cf. (I69l)). 

Fi = -2^2^' + G3Z + C4, Gi = 0, (73) 

where Co_4 and B12 are real constants. The conditions H = and A = exclude type 
III and N as special cases. The solution corresponding to ( 1721) with Ci = — | coincides 
with a special case of Kerr and Debney's type II solution [12] (p. 608), and hence is 



not new (see [Appendix D for more details). Also, by setting Co = 0, C2 = + C4 



and Bi = — IC1C3, one can obtain ( 173|l from (I72|) in the limit Ci — )■ 0, whereas the two 
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Killing vectors (151 p become commuting. For a generally non-vanishing A, we obtain 
from fH3l) an equation for J{z): 



J" 



= ^ _ 2AJJ" - yA(J')' - 2(AV2 + Cl)f 

+ 8{J'Y^\xp (^-3A j Jdz^ [5isin(3Ciz) - ^2 cos(3Ciz)], 



(74) 



which generalizes the equation for type II. By introducing J = K' (or — 3AJ = 
K' /K, etc.), one can transform this equation to a fourth-order ODE for Ki^z). However, 
we do not have any type II solution with A 7^ for this ODE. 

9. Conclusions 

With the real coordinates {x^z^u.r} and A(C) = 2 in the ansatz we present here 
our new class of vacuum twisting type II metrics admitting two Killing vectors: 

g = o^'im [dCdC + A (dr + WdC + WdC + HA)] , 

COS 2 / 

with C, = X + iz and 

e^i^dn - 2 [/exp (/ Fsd^) d^] dx 



A 



exp (/ F2(\.z) 



l(2F[ 
2 V ^1 



U 



^ ) - AF^ 



F2 + iCiJ (e-'^ + l), 

2(F3(cosr + 1) + F4sinr) 



(cosr -|- 1) AF? cosr, 

6 



(75) 



where Ci is an arbitrary real parameter and the real functions Fi_^{z) satisfy 
= -F'; + F,F[ + ^-KFl + i(3F^ - F| - Cl)F, + ^, 
= (i/' - 2F2iJ)' - 2F2(iJ' - 2F2iJ) + AClH, 
F3 = exp j F2dz^ [Bi sin(3Ci2) - B2 cos(3Ciz)], 

F4 = exp j F2dz^ [B2 sin(3Ci^) + Bi cos(3Ci^)], 

with the function H{z) defined by 

H = F'iFi - {F[f - AF^ - Fl^Fl 

One can determine Fi^2 from the first two equations of (175|) . which are in fact a pair 
of ODEs for Fi and K = j F2dz. This allows the metric to have at most eight real 
parameters including A and Ci (see discussion in Section 7), and its two Killing vectors 
are 



9,, 



- Ciudu- 
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Most of the previously known twisting vacuum solutions with two Killing vectors, as 
presented in Chapters 29 and 38 of [T^ , have been shown to belong to this class, the only 



exception being the general case of Kerr and Debney's type II solution (see Appendix D ). 
Additionally, for H = 0, the system ( 175!) can be reduced to a single ODE (1741) . or even 
be fully integrated when A = 0. This leads to the discovery of a limiting solution (l69ll 
(type II with two commuting Killing vectors) of a special case of Kerr and Debney's 
type II solution, which we believe has not been discussed or published before. Despite 
all these special solutions with maximally four parameters, the general solution of (175|) 
is still quite unknown. We believe that this problem poses a major challenge. 

Altogether, we hope that this work may provide a platform for all types of 
twisting algebraically special solutions to be studied in a connected and unified manner, 
given the history that many of those known solutions were derived by quite different 
approaches or special assumptions. For future research, this new class of metrics may be 
further examined to study important issues such as the cosmic no-hair conjecture, the 
asymptotic stability of the Kerr solution [M], and the formation of rotating black holes 
that might be described by certain solutions of ( 175|) (see, e.g., [351 ES] and references 
therein) . 

Appendix A. Type D solutions: \l/3 = \Ef4 = 

For simplicity, we consider type D solutions with \I'2 7^ and \l/3 = \E'4 = 0. In order to 
acquire ^3 = with ^^3 depending on r, we need, at least, for the coefficient of e^^'' to 
vanish in (ITTl) . 

= 2dlogp + c, (A.l) 
which, through the ansatz ( l42l) . adds two more equations to the system: 

Ci = 0, F[ = (A.2) 
Thus we can simplify the original 0431146 p to 

4 n 4Fc! 

= -|aF,^ - ^, (A.3) 

i^3 = 3F2F3, (A.4) 
Fi = 3F2F4, (A.5) 

with the last equation ( H6l) being automatically satisfied by the equations presented 
above. The system f lA.2tlA.5]) can be fully solved with the general solution 



Pi = C3 sec I — — — ) , ^2 = -pr tan 



C2 / G2 V C2 ^ . . . 

C|(3 + 2AC|C|) 6/z + Co\ ^ „ + ^ ^ 

^ ^ sec^ I — ^ I , ^4 = ^4 sec'' ' 



Remarkably, though we have only started with one extra condition ( lA.ll) . it is enough 
for the solution flA.6P to be of type D, i.e., that we have 



-C|(3 + 2AC|C|) + 6iC|C4, ,,3 T T n 
^2 = ^^^^^ ^(e" + l)^ ^^ = ^^ = 0. 
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With A{() = 2 and (lA.6p . the resulting metric can be written as 



CI 



° 2cOs2(f)cOs2(^) 

with ( = X + iz and 



[dCdC + A(dr + m)] 



A = cos 

n = 



2 / Z + Cq 



C2 
1 



du — 2C2tan(—^ — - 



dx 



cos 



C2 ' 



r \ / COS r 2C4 sin r \ 1 ^2 a / ^ \ 

cos- 1 2 j ( + ^4 1 + qC'|A(cos 2r + 2 cos r) 



'2 ^3 / 

One can immediately remove the parameter Cohj z + Cq z. Besides the cosmological 
constant A, the metric contains three real parameters C2, C3 and C4. 



Appendix B. Type D solutions from classical symmetries 

Here we list a number of special solutions one may encounter when searching for group- 
invariant solutions from classical symmetries [21] of the system fl43H46l) or its various 
special cases. Incidentally, all these solutions turn out to be of type D, even though 
the condition 2\E'2 = 3\E'3\E'4 is never used in their derivation. As expected for type 
D solutions, they all have Ci = 0, which means that the two Killing vectors ( ISTl) are 
commuting. In what follows, C2-4 are real constants. 

We start with solutions with F3 = 0. For A = 0, we have 



Fi 



C2Z 



F, = - 



F2 = 2C2, 
3 



0, 



F4 = C4e' 



6C22 



For A 



Fo = — , 



F3 
F3 = 

z 

and cf. dSl)) 

5 

S'z 
2 

z 

In case of F4 = 0, they all become type O solutions of fl54l) or fl63|) . 
For solutions that admit non-vanishing F3 and A, we have 

Fi = C2 7^0, F2 = 0, F3 = -iC|A, F4 = C3 



-s^ < 0, we have fH 
3sz^ 

^'-2;~z^ 



0, 



0, 



F^9l. 



Fi 



F, = -- 



3C^ + 2C^A 



C3 



z z bz° z° 

Note that the first solution above is given by constants, which is a consequence of the 
translational invariance {z — )■ z + Cq) of the system ( 143H46|) . Lastly for A = 0, we obtain 



Fi = C,e 



C2Z 



\C2 



1^4^2 4C2 2 
9<-^3'-^2^ , 



C4e 



4C2 2 
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Appendix C. Comparisons of type II solutions with Ci = 

Here we compare the three type II solutions from Case 2, i.e., (1671) . ( I68l) and ( l69l) . and 
show that fl69l) is different from the other two. Generally, to see that two twisting type 
II vacuum metrics are different, i.e., not being related by a coordinate transformation, 
it is sufficient to show that their CR structures along the shearfree null congruences are 
not equivalent [6]. This can be decided by evaluating the six Cartan invariants [H [T3], 
which are the same only for two equivalent CR structures. Among these invariants, the 
first one, in Cartan's original notation, is given by 



«(C,C) 



5rd(^r + rS^r + 8crr 
f ■ \/(rfy 



r = I (9^/ + 2c/) , / = —dc_d^c — cd^c, 

which only relies on the function c = c{(,() (hence F2{z) from the ansatz ( H2l) ). The 
use of a alone will be adequate for our comparison. 

For simplicity, we consider the special case of ( l69l) with C2 = = 0, i.e., 

2 2 

Fi = —2B2Z , F2 = -, F3 = — i?2-^i , F4 = BiF^ . 
This solution is still of type II but has a constant invariant a given by 

14 

Yet another case with a being constant is when i?2 = in fl69l) . in which case the 
solution is of type D and 



^2 = -I-VIO. 

5 

This same quantity a calculated from Lun's solution (j67j), however, is generally a 
function of z and only becomes a constant when two of the three free parameters £'1,2 
and M vanish, i.e., that we have 
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a^ = ll^ for E, = E2 = 0, 
10 12, 



a" 



72(4 - yi3) for El = M = 0, 



V2(4 + v^) for ^2 = M = 0, 



none of which is equal to those of ( l69l) . As for Demiahski's solution ( 168|) . one can see 
(using Maple) that its invariant a is never a constant (even when 6 = for type D; 
particularly for the NUT solution with 6 = a = 0, a is not defined due to f = 0, and the 
corresponding CR structure is hyperquadric [7]) within the full range of the parameters 
M, a and b. Therefore we conclude that the solution fl69|l is different from fl67|) and 
(EHD. 
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Appendix D. Kerr and Debney's type II solution 

The solution by Kerr and Debney [13] (see also [I2] P- 608) admits two non-commuting 
Killing vectors and reads 

P, = l, L = A,CV^- + A,CC^^ 

iMs -m, = 2Ai(l + a)C, A = 0, 

with Reo" = —3 and Ai^2 complex constants. The special case with a = —3 can be 
captured by the ansatz P2|) . and it corresponds to (-Ei_4 real) 

^(C) = C, C=IC|e^^^ ^1 = -^ 
Fi = E3 sin(f ) - cos(f ) + 2Ei sin(f ) - 2E2 cos(f ), 
F2 = F[/F,, F-s = -F^ (El sin(f ) - E^ cos(f )) , 
F4 = Fi'(E2sin(f) + Ei cos(f)), 

Ai= Ei + 1E2, ^2=^3 + 1^4- 

Modulo some redefinition of parameters, this is the same solution as ( 1721) and (1711) with 
Ci = — |. It is not clear how to make such a conversion for the general case of fID.ll) . 
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